Three-dimensional skyrmion states in thin films of cubic helimagnets 
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A direct three-dimensional minimization of the standard energy functional shows that in thin films 
of cubic helimagnets chiral skyrmions are modulated along three spatial directions. The structure of 
such 3D skyrmions can be thought of as a superposition of conical modulations along the skyrmion 
axis and double-twist rotation in the perpendicular plane. Numerical solutions for chiral modulations 
demonstrate that 3D skyrmion lattices and helicoids are thermodynamically stable in a broad range 
of applied magnetic fields. Our results disclose a basic physical mechanism underlying the formation 
of skyrmion states recently observed in nanolayers of cubic helimagnets. 



The Dzyaloshinskii-Moriya (DM) interactions arising 
in noncentrosymmetric magnets as a result of their 
crystallographic handedness, induce long-range modu- 
lations with a fixed sense of the magnetization rota- 
tion [1]. It was also established that DM interac- 
tions provide a unique mechanism to stabilize two- and 
three-dimensional (2D and 3D) localized states {chi- 
ral skyrmions and hopfions) and their bound states 
(skyrmion lattices) [2-4]. 
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Figure 1: (color online). Magnetic configurations in 
basic chiral modulations arising in cubic helimagnets: 
(a) an isolated skyrmion with double-twist modulations 
perpendicular to the applied field, (b) transversally dis- 
torted helices {helicoids), (c) a hexagonal lattices of chiral 
skyrmions (a), (d) longitudinal distorted helices {cones). 



Skyrmion lattices and isolated skyrmions have been re- 
cently discovered in nanolayers of cubic helimagnets [5- 
8] and in monolayers of magnetic metals with induced 
DM interactions [9]. Chiral skyrmions are highly mo- 
bile nanoscale spots of reverse magnetization that can be 
utilized in ultra dense recording technologies and in the 
emerging spin electronics [6, 10]. Importantly, in the ma- 
jority of nonlinear field models static multidimensional 
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solitons are unstable and collapse spontaneously into 
topological singularities [11]. This singles out condensed- 
matter systems with broken chiral symmetry as a special 
class of materials, where chiral skyrmions can be formed 
as static, mesoscopic solitons and can be manipulated in 
broad ranges of the thermodynamical parameters [2, 3]. 

The experimental observations of skyrmions and other 
chiral modulations [5, 7, 8, 12] are in an excellent agree- 
ment with results of the phenomenological theory based 
on the classical Dzyaloshinskii model [1, 3, 12, 13]. Par- 
ticularly, theoretical analysis shows that induced uniaxial 
anisotropy stabilizes skyrmions and helicoids in thin lay- 
ers of cubic helimagnets [12, 14]. In bulk cubic helimag- 
nets, these chiral modulations exist as met ast able states 
[3, 14]. Recent experiments in MnSi and FeGe confirm 
that they are stabilize by an induced uniaxial anisotropy 
[7, 8, 12]. On the other hand, nothing is known about 
a role of induced magnetic anisotropy in mechanically 
thinned nanolayers of cubic helimagnets where direct ob- 
servations of skyrionic states have been reported [5, 6]. 
Generally, the exact magnetic structures of skyrmions 
and other modulated states arising in confined cubic he- 
limagnets are still not resolved and physical mechanisms 
underlying their formation and stability are not com- 
pletely understood. 

In this paper we provide the first three-dimensional 
calculations of chiral modulated states in thin layers of 
cubic helimagnet films. We show that below a critical 
thickness, a double- twist skyrmion gains additional mod- 
ulations along its axis. The calculated phase diagram for 
3D skyrmions and competing phases indicates a broad 
range of the films thickness, where skyrmionic states are 
thermodynamically stable. 

Within the phenomenological theory introduced by 
Dzyaloshinskii [1], the magnetic energy density of a cu- 
bic non-centrosymmetric ferromagnet can be written as 
[1, 15] 

w = A (grad m)^ + i:> m • rot m - H • m M. (1) 

m = (sin^cos?/;; sin^sin?/;; cos^) is the unity vector 
along the magnetization M and H is the applied field. 
Energy density functional (1) includes only the principal 
interactions essential to stabilize skyrmions and helicoids: 
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the exchange stiffness with constant A, DM couphng en- 
ergy with constant and the Zeeman energy. We ne- 
glect less important energy contributions such as intrin- 
sic and induced magnetic anisotropics or stray-fields. In 
contrast to common magnetically soft materials where 
magneto-dipole interactions play an important role [16], 
stray-field effects in noncentrosymmetric magnets are 
found to be weak due to the stabilizing influence of the 
DM interactions [10]. 

DM interactions favour spatial modulations of M with 
a fixed rotation sense [1]. In an infinite sample, isolated 
skyrmions (Fig. 1 a) are described by solutions of type 
[2, 3], 

e = e{p), v^ = (^ + ^/2. (2) 

Here we write the spatial variables in the terms of 
cylindrical coordinates r = (pcos(/?; psin(/?; 2;). Eq. 
(2) describes 2D skyrmions as axisymmetric tubes with 
"double-twist" modulations in the (x, y) plane and a ho- 
mogeneous distribution along the skyrmion axis (z-axis) 
(Fig. 1 a) [3, 17]. The solutions for one-dimensional 
modulations include distorted helices with the propaga- 
tion direction perpendicular to the applied field, helicoids 
(Fig. 1 b), and longitudinal distorted helices modulated 
along the field, cones (Fig. 1 d). The analytical solu- 
tions for helicoids describe a gradual unwinding of chiral 
coils in an increasing magnetic field [1]. The equilibrium 
parameters for the cone phase are given by [15] 

cos6> = IHI/Hd, ip{z) = V^cone(^) = 2itz/Ld, (3) 

where the period of the helix, and the saturation field, 
Hd are expressed as 

Ld = ^^A/\Dl Hn = Dy{2AM). (4) 

We introduce the reduced energy density of a modulate 
phase as e = (w) /[D'^ /{AA)] where (w) is the energy den- 
sity averaged over the modulation period. Then the equi- 
librium energy density of the cone phase can be written 
as 

econe = - [1 + [H/Hof] (5) 

In an infinite sample, the cone (3) corresponds to the 
global minimum of model (1) in the entire magnetic field 
range below the saturation field (4). Helicoids and 
skyrmion lattices only exist as metastable states [14]. 

Contrary to lower symmetry noncentrosymmetric fer- 
romagnets [2] , in cubic helimagnets DM interactions pro- 
vide three equivalent modulation directions. However, 
2D chiral skyrmions investigated in [3, 13] have only 
two (in-plane) modulation directions while along the 
skyrmion axis the structure remains homogeneous. Thus, 
the question arises: is it possible to lower the skyrmion 
energy by imposing additional chiral modulations along 
the skyrmion axis? Our calculations based on rigorous 
solutions of micromagnetic equations show that the for- 
mation of 3D skyrmions actually occurs in the thickness 
range below a critical thickness. 



As a model we consider a film of thickness L, infinite 
in X and y directions and bounded by parallel planes at 
z = ±L/2. Numerical solutions for chiral modulations 
are derived by direct minimization of functional (1) with 
free boundary conditions along the z-axis and periodic 
boundary conditions in the (x, y) plane. This yields the 
equilibrium spatially inhomogeneous distributions of the 
magnetization vector m in the layer as functions of the 
three spatial variables and the two control parameters, 
reduced magnetic field, H/Hjj and confinement ratio^ 
L/Ld. 

A three-dimensional minimization of nonlinear sigma 
models needs huge computing efforts and is usually exe- 
cuted on high-performance supercomputers [18]. Modern 
graphics cards offer a promising alternative. However, 
this method requires special hardware-oriented highly 
parallel algorithms [19]. We have developed a special 
algorithm for CUDA architecture for NVIDIA graph- 
ics cards. This architecture has been found to be effi- 
cient to solve complex micromagnetic problems [20]. We 
use a nonlinear conjugate gradient method to minimize 
energy functional (1). The results of the minimization 
have been checked by the compatibility with the Euler- 
Lagrange equations. The results of numerical minimiza- 
tion of model (1) are presented in Figs. 2, 3, 4 and as 
video [21]. 




Figure 2: (color online). 3D skyrmion lattice: (a) Distri- 
bution of the magnetization in the skyrmion core demon- 
strates chiral conical modulations along the cell axis. For 
clarity, the sizes along the z-axis are magnified; (b) Cal- 
culated distribution of the magnetization in the skyrmion 
cell for top and bottom layers in a film with thickness 
L/Ld = 0.25 and in the applied field H/Hd = 0.2. 

Isolated and embedded skyrmions. The equilibrium 
magnetic configurations in a 3D skyrmion lattice strongly 
differ from axisymmeteric 2D skyrmions arising in infinite 
helimagnets (2). The magnetization vector m depends on 
all three spatial components, (p, z), and exhibits com- 
plex three-dimensional modulations ((Fig. 2) and [21]). 
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It was, however, found that a simphfied ansatz 

e = e{p), ^(0,^) = ^ + 7r/2 + ^?(z) (6) 

with iIj{z) = V^cone (3) provides a nice approximation of 
the solutions for isolated and bound 3D skyrmions. 
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Figure 3: (color online). Calculated equilibrium ener- 
gies of isolated skyrmions with ansatz (6) (solid line) 
and for those that are homogeneous along the z-axis (2D 
skyrmions, dashed line) as functions of the reduced field 
H/Hd in a film with reduced thickness L/Ljj = 0.25. 
Inset shows the equilibrium energies of the skyrmions as 
functions of the reduced thickness at a fixed applied field 
H/Hd = 0.4. 

Eq. (6) describes double-twist rotations in the (x^y) 
plane and conical modulations (3) along the skyrmion 
axis (Fig. 2 a). Remarkably, the period of conical mod- 
ulations in the skyrmion coincides with the helix period. 
We introduce the reduced equilibrium energy of an 
isolated skyrmion as 



{w - wo)pdp, (7) 



where wq = —HM is the energy density of the saturated 
state. In Fig. 3 The energies Ws for isolated skyrmions 
with ansatz (6) and for solutions homogeneous along the 
skyrmion axis (2D skyrmions) are plotted as functions of 
the applied field in Fig. 3 (main figure) and the reduced 
film thickness (inset). These results show that in a broad 
range of the control parameters H/Hjj and L/Ld-, the 
modulation along the skyrmion axis reduces its energy. 
3D modulated skyrmions should not be confused with 
three-dimensional topological chiral solitons (hopfions) 
[4]. 

Helicoids and cones. In thin films, helicoids become 
inhomogeneous along the film thickness [21]. However, 
we failed to find a simple ansatz to describe these dis- 
tortions. Conical states propagating perpendicular to 
the layer surfaces are compatible with the film geometry, 
and the free boundary conditions impose no constrictions 
on these modulations. As a result, the solutions for the 



cone phase in the film and the equilibrium energy den- 
sity equal to those in a bulk sample (3), (5) [14, 15]. 
Conical states (3) can exist even in films with thickness 
smaller than the helix period { L < Ljj) (4). In this case 
AV^ = |V^(L/2) - ^|J{-L/2)\ = 2tiL/Ld < 27r. 

The calculated equilibrium energy densities of heli- 
coids (eh) and 3D skyrmion lattices (cgk) have been com- 
pared with the cone energy density (econe)- The re- 
sults are plotted as functions Aeh = (eh — econe) and 
A^sk = (csk — Gcone) vcrsus the applied field together 
with the corresponding results for ID helicoids and 2D 
skyrmion lattices arising in bulk cubic helimagnets [1, 3] 
(Fig. 4 a). In the latter case the cone has the lowest 
energy in the entire region below the saturation field, 
Hd (3). This situation retains in sufficiently thick films 
(see Fig. 3 Inset). However, the energy balance between 
the chiral phases drastically changes in thin layers below 
the critical thickness. Here the formation of 3D modula- 
tions is resulted in a large reduction of the helicoids and 
skyrmion lattices and leads to their energetic advantage 
over the cone phase. (Fig. 4 (a)). It means that in suffi- 
ciently thin cubic helimagnets fflms the onset of conical 
modulations in helicoids and skyrmion lattices provides a 
speciffc mechanism to stabilize these chiral modulations 
in a broad range of applied magnetic ffelds (Fig. 4). 
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Figure 4: (color online), (a) The differences Ae be- 
tween the energies of skyrmion lattice and the cone phase 
(solid, blue) and between the helicoid and the cone phase 
(dashed, green) as functions of the applied field for 2D 
and 3D modulations. Functions Ae{H/HD) are plotted 
for two values of reduced film thicknesses L/Ljj = 0.25 
and 0.5 in comparison with results derived for ID heli- 
coids [1] and hexagonal skyrmion lattices homogeneous 
along the skyrmion axis (2D skyrmions investigated in 
[3]). (b) Magnetization curves for 3D modulated struc- 
tures in a film with reduced thickness L/Ljj = 0.25 (solid 
lines) in comparison with the magnetization curves for 
ID helicoids and 2D skyrmion lattices [1, 3]. 



An elementary analysis of modulations (6) allows to 
understand basic physical mechanisms underlying the 
formation of 3D skyrmions. 
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Inserting ansatz (6) into (1) yields the following ex- 
pression 



w = A e 



sm^O\ ^ T/ N sin 2(9 

P ■ ) + p cos [ 0, + ^ 



H McosO^ sin^ (9 



(8) 



In functional (8) three first terms are responsible for 
the formation of double-twist modulations in the (x, y) 
plane [3] and Wz{z) describes modulations along the 
skyrmion axis. With conical modes = '^cone (3), 

the function D{z) and energy density Wz can be written 



D{z) = Dcos{27rz/LD) , = - sin^ / {4A) . (9) 

Energy density (8) includes two different DM terms. 
The first, with factor D{z), favours double- twist rota- 
tions in the (x, y) plane, the other induces cone-mode 
modulations along the z— axis. The energy gain from 
the double- twist modulations oc D'^{z) reaches the max- 
imum at tp{z) = 0. The cone modulations increase the 
double-twist energy. On the other hand, these yield a 
negative energy contribution Wz (9) into the total en- 
ergy. The equilibrium 3D skyrmion patterns are formed 
as a result of the competition between different DM in- 
teraction terms. 

In the calculated 3D skyrmion textures angle i/j equals 
zero in the layer center and rotates with a same rotation 
sense to the maximal value at the surfaces, |?/;(±L/2)| = 
7tL/L]j. Importantly, for ip > 7r/2 the factor D{z) be- 
comes negative. It means that cubic helimagnets films 
with thickness larger than Ld/2 include regions where 
function D{z) becomes negative corresponding to the 
double-twist modes with energetically unfavorable sense 
of rotation. This explains the existence of the critical 
thickness for 3D modulations and allows to estimate its 
value. As a rule of a thumb, the critical thickness can be 



estimated as a half of the helix period Ld (see e.g. Fig. 
3 Inset). 

The first direct observations of chiral skyrmions have 
been reported in a mechanically thinned cubic helimag- 
net (Fe,Co)Si with thickness L = 20 nm [5]. In this 
material Ljj = 90 nm, and, thus, confinement ratio 
L/Ld = 0.22. This is quite below the critical thickness, 
and ,thus, 3D modulations are expected to exist in this 
film. In Fig 4 (b) we present the calculated magnetization 
curve for a film with confinement ratio L/ = 0.25 close 
to that of the film investigated by Yu et al. [5]. This in- 
dicates the stability regions for helicoids and skyrmions, 
the first order transitions between helicoid and skyrmion 
lattices, and the second order transition of the skyrmion 
lattice into the saturated state. 

Earlier an alternative mechanism to stabilize chiral 
skyrmions in cubic helimagnets has been proposed by 
Butenko et al. [14]. According to this paper the dis- 
tortions of a cone phase imposed by uniaxial magnetic 
anisotropy strongly increase its energy. As a result, 
chiral skyrmion lattices and helicoid become thermody- 
namically stable in a broad range of applied magnetic 
fields (Fig. 3 in [14]). In our recent studies, the cone 
modes remain undistorted. The energetic advantage of 
3D skyrmion lattices and helicoids is gained because 3D 
modulations grant a larger reduction of the DM energy 
than single-direction modulations in cone modes. New 
experiments are needed for detailed investigations of in- 
duced magnetic anisotropy and to resolve the exact mag- 
netic structures in confined cubic helimagnets. These will 
allow to understand a role of unaxial distortions and 3D 
modulations in the formation of skyrmion states in these 
systems. 

In conclusion, numerically calculated 3D solutions for 
chiral modulations in nanolayers of cubic heliimagnets 
show that in the films below the critical thickness (es- 
timated as a half of the helix period) chiral skyrmions 
and helicoids are sufficiently inhomogeneous along three 
spatial directions. Such 3D modulated skyrmions and he- 
licoids modulated along the layer thickness are thermo- 
dynamically stable in a broad range of applied magnetic 
fields. 

We thank T. Monchesky for critically reading the 
manuscript and fruitful advice. 
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